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g^ ■ Abstract 



^ . We estimate freezeout conditions for s, c, and b quarks in high energy nuclear 

collisions. Freezeout is due either to loss of thermal contact, or to particles 
"wandering" out of the region of hot matter. We then develop a thermal 
recombination model in which both single-particle (quark and antiquark) and 
two-particle (quark-antiquark) densities are conserved. Conservation of two- 
particle densities is necessary because quarks and antiquarks are always pro- 
duced in coincidence, so that the local two-particle density can be much larger 



O ! than the product of the single-particle densities. We use the freezeout condi- 

^^ ' tions and recombination model to discuss heavy resonance production at zero 

baryon density in high energy nuclear collisions. 
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I. INTRODUCTION 

One of the striking features of high energy proton-proton, proton-nucleus and nucleus- 
nucleus collisions is the production of heavy resonances. Since these resonances appear as 
a consequence of raising the collision energy, it is natural to ask, "Can heavy resonances 
be used experimentally to signal a deconfined phase of matter?" Enhanced strangeness 
production was in fact proposed as one of the first potential signals of the Quark Gluon 
Plasma (QGP). |I|] i/ip production can be looked upon as changing the focus to the 
heavier charm (c) quark from the much lighter strange (s) quark. Here we aim at a unified 
treatment of heavy quark composites, extending the scope to include also the bottom (6) 
quark. 

In view of the very different masses of the s, c, and h quarks, and of the mesons that 
contain these quarks, it is reasonable to assume that the various flavors will stop interacting 
with the expanding collision system at different times and thus carry "frozen out" infor- 
mation about different stages of the time evolution, reflecting their last interactions. This 
calls for a detailed examination of the freezeout conditions for each flavor. The concept of 
different freezeout times for different species is well known in studies of the Early Universe, 
and has also been discussed in the context of nuclear collisions at lower energies. 

Our physical picture assumes thermal equilibrium, and chemical equilibrium between 
resonances with the same heavy quark content (but not between heavy quarks and anti- 
quarks) at high temperature. We allow resonances to depart from chemical equilibrium 
after freezeout. In determining the freezeout chemical potentials we pay particular atten- 
tion to the local quark-antiquark {QQ) pair density. Since quarks and antiquarks are always 
created in pairs, the local pair density may be signiflcantly higher than the product of the 
single-particle densities. We take this important effect into account in terms of a pair chem- 
ical potential (see Section |T|). This pair chemical potential is the main difference between 
our recombination model and others that have been used recently to study strange [^§ and 
heavier exotic resonances. 

In Section we discuss two conditions which may signal that heavy resonances cease to 
interact with the hot matter, and estimate freezeout temperatures and proper times for s, 
c, and h quark resonances in selected AA collisions at CERN's Super Proton Synchrotron 
(SPS, for which ^/s = 20 GeV/nucleon) and future Large Hadron Collider (LHC, i/i = 7 
TeV/nucleon), and at Brookhaven's future Relativistic Heavy Ion Collider (RHIC, y/s = 200 
GeV/nucleon). We obtain a rich variety of freezeout scenarios, depending on the properties 
of the transition from the hadronic phase to the high temperature (deconflned) phase. We 
also calculate the mean distance a heavy quark moves before freezeout and the mean number 
of interactions it has; the mean number of interactions is high enough that thermal recombi- 
nation models can be expected to work for heavy ion collisions at or above SPS energy, but 
not necessarily for collisions of lighter ions. In Section |TI| we describe our thermal recom- 
bination model. We then give an illustrative application to J ftp suppression, and calculate 
freezeout densities and the fractions of c and b quarks that end up in cc and bb resonances. 



Finally, in Section IV we summarize our work. 



II. FREEZEOUT CONDITIONS 

In this section, we estimate the freezeout conditions for the different heavy quark flavors. 



We consider two mechanisms for freezeout. In subsection |IIA| we discuss freezeout due to 
loss of thermal contact, which we refer to as "contact" freezeout. In subsection [11 B| we 
discuss freezeout by emission, which we estimate by treating the quarks as random walkers, 
and hence refer to as "wandering" freezeout. We then summarize the results of the two 
freezeout scenarios in subsection [II Q 

The s quark mass, rus, is approximately equal to the quantum chromodynamics (QCD) 
transition temperature, Tc ~ 150 — 200 MeV, while the c and b quark masses, respectively 
rric and rrih, are much greater. [We use rus = 0.2 GeV, rric = 1.5 GeV, and mi, = 5 GeV.] 
Thus, freezeout conditions for c and b quarks are slightly different than those for s quarks. 
Because of this, we estimate freezeout conditions under each scenario for c and b quarks 
together, and then for s quarks. 

A. Loss of thermal contact 

One possibility is that freezeout occurs when the thermal contact between the heavy 
quark and the hot matter ends. In this case, the temperature at freezeout, Tf, is determined 
by setting the mean free path equal to the mean distance to leave the hot matter. If the 
hot matter forms an infinite cylinder of radius Rf at freezeout, the initial location is equally 
likely to be anywhere within the cylinder, and any initial direction of motion is equally likely, 
the mean distance to leave the hot matter is 4Rf/3. The thermal heavy quark mean free 
path is 

where vq is the velocity and Tq is the collision rate. For heavy quarks in deconfined matter, 
we take 
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^Q (T) = ^^jVmQT, (4) 

where we use Tq from Ref. 0] with A^ = 3 colors, Casimir factor Cf = 4/3, Nf = 3 light 
flavors (m, d, s) and strong coupling constant g = 2. We have dropped terms of order 
Vq ln(l/fQ) in these equations, as these are small and do not affect our results much. 
In the coexistence region, the mean free path is 
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where x is the fraction of matter in the deconfined phase. We estimate mean free paths 
in hadronic matter at T^ by assuming that the density of scatterers is proportional to the 
entropy density in the two phases, and that cross sections remain constant: 

AW(T,) = ^A('^)(T,), (6) 

where z/ = Sd{Tc)/sh{Tc) is the ratio of entropy densities in the two phases at Tc. This is 
equivalent to assuming that the number of scatterers is proportional to the entropy, and that 
the phase transition dilutes the scatterers without affecting their cross sections; the extra 
factor of 1/2 occurs because quarks must travel in mesons (which we treat as two quarks) 
in the hadronic phase. We thus obtain 

Ag) = -i^ _ Ag) (Te) . (7) 

For hadronic gas, we expect the mean free path to be 

Aq (T) - -^, (8) 

where p is the density of light particles (with m <ti T) and aq is the mean cross section. The 
thermal velocity is approximately the same for heavy quarks and mesons containing heavy 
quarks, so at Tc we find 

assuming continuity of Xq, and that p is given by the density of a massless pion gas, 

P = ^T^- (10) 

For Tc = 150 — 200 MeV and u = 5 — 10, aq = 2 — 8 mb, in reasonable agreement with 
currently used values for heavy meson cross-sections. However, different heavy quark mesons 
may have very different mean free paths, so it is not clear that all meson species containing 
any given heavy quark flavor will freeze out at the same time. This issue is beyond the scope 
of this paper, so we leave it for future work. 

We extrapolate to lower temperatures by assuming that ctq is constant and p oc T^, so 

^^q\t) = ^X^\Tc) [^) . (11) 

This estimate of Xq is reasonably accurate near Tc. However, for T <C m.,^ the mean free 
path will grow exponentially with T~^ as the density drops, so freezeout will occur before 
T drops much below rrij^. 

For freezeout occurring from deconfined matter, we obtain 

^Rf = ^-^PI^QTf ^Tf = ^^ruQ^Rj' ^ Qm^' Rj'. (12) 



Freezeout occurs for T > Tc if 

Rf < ^(mQ^.)-/^^3(mg^,)-/^ (13) 

"Contact" freezeout occurs from the coexistence region as long as 

i?/ < ^ ^ imQ^r'^' ^ 1.5 u (rngT,)-^/^ , (14) 

and occurs from normal hadronic matter for larger values of Rf. Note that in all equations 
above Rf can depend on Tf or Tf, so that the equations are valid even with transverse 
expansion of the hot matter as long as Rf is treated properly. 

The s quark mass is small, so its velocity can be taken to be unity, and the mean free 
path is 

\if)^T) = —^ = 2T-\ (15) 

where we obtain r^"^-* for light quarks from Ref. 0. We then find that "contact" freezeout 
of s quarks from deconfined matter occurs when 
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Tf = -RJ\ (16) 



2 
Freezeout thus occurs at or above T^ if 

Rf < It-\ (17) 

We estimate the mean free path for an s quark in the coexistence region just as we did for 
the c and b quarks. The thermal velocity in hadronic matter is still approximately unity near 
Tc, as the mass of the lightest meson containing an s quark is not too much greater than the 
transition temperature. If we again assume that the number of scatterers is proportional to 
the entropy, and cross sections are constant, we obtain Eq. (|^) as before. This is equivalent 
to assuming a cross section 

a, = -^^ T-^ ~ 3 - 10 mb, (18) 

which again is not an unreasonable value. We find that "contact" freezeout of s quarks in 
the coexistence region occurs if 

Rf<l^T^~\ (19) 

For larger values of Rf, the s quarks will freeze out from normal hadronic matter. 



B. Emission ("wandering") 

A second possibility is that freezeout occurs when the quarks leave the hot matter. This 
happens when the path length becomes equal to the mean distance to leave the cylinder. We 
calculate the path length by treating the quarks as random walkers, with thermal velocities 
and mean free paths. The rms distance moved by a quark (or antiquark) of flavor Q, tq, 
from creation at proper time tq ~ to freezeout at proper time Tf is given by 
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where vq and Xq are respectively the thermal velocity and mean free path for flavor Q, 
evaluated at T{t). We relate T and r by treating the deconfined matter as a non-interacting 
gas with 16 boson and 24 fermion degrees of freedom, assuming that entropy is conserved, 
M and that the final state can be treated as a massless Bose gas, so that 
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747]- 
ttRJtx ——T^ = 3.6dN/dy, (21) 

where dN/dy is the total rapidity density (charged plus neutral). As we neglect transverse 
expansion, we set Rf equal to the nuclear radius, 

Rf = 1.2 A^^^ im. (22) 

The assumption that Tq = is good for c and b quarks, which are only created early in the 
collision, but not necessarily for s quarks which are created as long as there is any deconfined 
matter. [Q However, we make this approximation for simplicity. 

"Wandering" freezeout of heavy quarks occurs in deconfined matter when 

^Rf) ^^^^^ dN/dy > ^^3 mg T^ i?J ^ 3.6 rug T<? RJ, (23) 

where Tc is the proper time at which the deconfined matter reaches Tc. Including transverse 
expansion increases Rf and decreases T for fixed r, making it less likely that c or b quarks 
freeze out before the hot matter reaches T^. The condition for "wandering" freezeout in the 
coexistence region is 



23687r2 
6561 {1 + (z^ - l)[(z/-l) - (z^ - 2) ln2]} 



dN/dy > .,.. n , ... .^... .^ .. o^ ,^ on ^Q ^' R% (24) 



if there is no transverse expansion. For smaller values of dN/dy, "wandering" freezeout 
occurs during the hadronic phase. 

For "wandering" of s quarks, we again take Eq. (120|), but with vq = 1. Using Xs{T) from 
Eq. (pISD, we find that s quarks freeze out from deconfined matter if 
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dN/dy > 1^ T^ R) ^ 16 T^ RJ. (25) 



"Wandering" freezeout from the coexistence region will occur if 
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dN/dy > /' , — -- T^ R) . (26) 

' ^ - 729{3 + 4(z/-l)[(z/-l)-(z/-2)ln2]} '^ ^ ^ ' 

For smaller values of dN/dy, "wandering" freezeout takes place from hadronic matter. Again, 
including transverse expansion would shorten the duration of the hot matter at any entropy 
density, so the above inequalities give solid lower limits on the values of dN/dy needed for 
"wandering" freezeout. 

C. Freezeout summary 

We have evaluated freezeout conditions for s, c, and b quarks in proton collisions at SPS, 
RHIC, and LHC energies. We find that the mean number of collisions before freezeout. 
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rdrTlTir)], (27) 
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is much less than unity for proton collisions at any feasible energy. Thus, the behavior 
of heavy quarks in proton collisions probably cannot be described by the quasi-equilibrium 
models used here, and so predictions for proton collisions are beyond the scope of this paper. 

In Table I, we show freezeout conditions for the various heavy quark flavors in selected 
AA collisions at SPS, RHIC, and LHC. We use four scenarios, with z/ = 5 and 10, and 
Tc = 150 and 200 MeV. This gives a reasonable range of values for Tc, while i/ = 10 is near 
the upper limit for a strong first-order QCD phase transition, and z/ = 5 gives a moderately 
strong phase transition. These values are compatible with lattice gauge theory results. 

The number of collisions per quark is only weakly dependent on u, T^, and flavor. For 
S+S collisions at SPS, Uf ~ 1, so quasi-equilibrium models are not ridiculous but are also 
not very well justified. For Au+Au collisions, Uf is significantly larger than unity, and 
thus quasi-equilibrium models are at least self-consistent even at SPS. This is in general 
agreement with Parton Cascade Model (PCM) results for Au+Au collisions at RHIC and 



LHC, [|IOl which show that all quarks follow approximately hydrodynamic behavior, with 
mean momenta following the temperature of the surrounding matter. 

None of the quark flavors freeze out above T^. The s quarks freeze out before the end 
of coexistence in any collision if i/ = 10. For z/ = 5, they may freeze out in the hadronic 
phase, especially if Tc is relatively high. In general, freezeout of s quarks typically occurs at 

The c quarks also freeze out in the coexistence region if z/ = 10, except in heavy ion 
collisions at SPS and RHIC if T^ is large. If z/ = 5, they freeze out below T^. Since these are 
upper limits, it is likely that c quarks also freeze out in the coexistence region if there is a 
strong QCD phase transition, but they freeze out from hadronic matter if the transition is 
weak. 

The b quarks may freeze out in the coexistence region in Au+Au collisions at LHC, but 
only if the transition is strong and T^ is relatively low. In all other cases considered, they 
freeze out far below Tc. Thus, accurate predictions of b quark physics in Au+Au collisions 
below LHC energy (and possibly at LHC, if Tc ~ 200 MeV) probably requires accurate cross 
sections for the various b mesons in hadronic matter, and possibly the use of a microscopic 
hadronic cascade code. 



The energy required for freezeout in the coexistence region will be lower in pA and 
AB collisions than in AA collisions. However, including transverse expansion will raise the 
energy thresholds, as might more realistic estimates of heavy resonance cross sections in 
hadronic matter, so the values given here are probably lower bounds for freezeout before the 
end of coexistence. 



III. THE THERMAL RECOMBINATION MODEL 

We use a thermal recombination model that is similar to the models of Refs. |§|-§]. The 
major difference is that we introduce a pair chemical potential to account for the fact that 
the quarks and antiquarks are always produced in pairs, so the local density of pairs may 
be much larger than the product of the quark and antiquark densities: 

pgl(x,x) » PQix)pQix). (28) 

We assume that the heavy quark flavors (s, c, and b) are far from chemical equilibrium before 
they stop interacting thermally, so that the densities of heavy quarks and antiquarks at 
thermal freezeout are separately conserved. The chemical potential for hadron species n is 
determined by the quark, antiquark, and pair chemical potentials, respectively /ig, /Iq, and 

/in = E {kPf^Q + kpfiQ + kl^^^f^Qo) . (29) 

Q 

Here /c^*^) (^i*^'*) is the number of quarks (antiquarks) of flavor Q, and k^^'^^ is the number 
of pairs (the smaller of k^^^ and kj^^). We assume that the number of hadrons of species n 
produced is proportional to its density at temperature Tf and chemical potential /i„. 

The chemical potentials are determined by conservation of quark and antiquark densities, 
respectively pg and pg, and the local two-particle quark-antiquark density, p^^. [We use 

the more compact notation p -^ = p Ux, x) for the remainder of this paper.] Since there is 
no chemical equilibrium between quarks and antiquarks, the constraint pq = —pq does not 
apply. The equations are 

Qn.kl^J r ^ p' dp ^ ^n^fcjy r P^ dp 

27r2 Jo e^E„^^P^-^^„^]/Tf _ i + f- 2tt^ 

_ ^ 9n,ki^J [ ^ p^ dp ^ ^n/^iy / - P- ap , V 

^ tf 27r2 io e[^"'>(P)-'^"f>l/^/ - 1 t;' 27r2 io e'^-z^^^"'^-/!/^^ + l' 
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^"6 4? f°° P^dp gnfk\^^ f°° p^dp 



"b 
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Pq tr 2772 7o g[s„jp)-M;j/T, _ ^ ' z_. 27r2 7o e'''"/(^)~'^''vl/^^ + l' 

fJ''n = f^n + Pqq, (33) 
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EAp) = ip' + mf'" 



(34) 



Here the sums are over bosons (n;,) and fermions (?^/), while m„ and (?„ are respectively the 
mass and degeneracy for hadron n. 

We estimate p^^ by again treating the quarks as random walkers. We assume that the 

Q and Q are distributed randomly over a sphere centered on the point at which the pair 
was formed, such that the rms distances from the center are tq. We thus obtain 



P, 



(2) 



PQ 
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PQ + 
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(35) 



where pg, fyQ and tq are determined at freezeout. We obtain the rhs of Eq. ( |32D by averaging 
the QQ pairs in hadrons over the same sphere to get the mean value of p j^ for the hadrons. 
For copiously produced quarks, p ^^ ~ PqPq^ so the third constraint equation (|34 ) just gives 



0. 



A. J/V' suppression 



One possible application of this model is to understand the phenomenon of J/ip suppres- 
sion. The prediction is that cc pairs will tend to form J/ip mesons less often in collisions 
in which QGP is formed, as a result of interactions of the pair with the deconfined matter. 
Thus, if QGP is formed in nuclear collisions, production of J/ip mesons will be suppressed 
relative to production of cc pairs. 

In this section, we show how suppression of heavy QQ mesons arises in our model. 
For illustrative purposes, we consider a case for which there are only three hadron species. 
Species 1 hadrons contain a single Q and have degeneracy gi and mass mi. Species 2 hadrons 
are the antiparticles of species 1, so they contain a single Q and have degeneracy gi and 
mass mi. Species 3 hadrons contain a QQ pair, and have degeneracy g2 and mass m,2 — 2mi. 
This is similar to the situation for c quarks, where the D and D could be regarded as species 
1 and 2 and the J/?/' as species 3. 

Because heavy quarks freeze out at Tf <^ tuq, we use Maxwell-Boltzmann statistics for 
the heavy hadrons. We take pq = pQ, because heavy quarks are always produced in pairs 
so the net density is always zero. We then solve the constraint equations, obtaining 



P3 

Pi 



1 + 
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gl 



^i2mi-m2)/Tf 
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PQ 



1/2' 



(36) 



P^ ) that the relative densities of the resonances depends only on Tf 
\l PQ- This remains true for more complicated resonance spectra 



It is clear from eq. 
and the local density, p|^ 
as long as all resonances have at most one quark and one antiquark, if the flavor is heavy 
enough to justify the use of Maxwell-Boltzmann statistics. These conditions are met by the 
flavors c and 6, so studies of the relative production of c and h resonances (except for exotic 



9 



multi-heavy hadrons 0) will be insensitive to pc and pb- Thus, the only way to determine 
pc and Pb (in the absence of exotic multi-heavy hadrons) is through comparisons of heavy 
and light resonance production. Since there are many multi-strange baryons, ps can be 
determined solely by studying strange resonances, as long as these multi-strange baryons 
are detected. 



B. Preezeout density estimates 

We estimate freezeout densities from rapidity densities, staying as close as possible to 
predictions supported by data. In the absence of transverse expansion, 

P« = Pq = 7^J7^' (37) 

and we obtain p j- from eq. (^). The predicted freezeout densities are shown in Table II; 
for the interested reader, we explain our estimates of dNq/dy below. The main point of 

(2) 

interest is that p J^ ^ PqPq fo^ ^ quarks in all colhsions considered, and for c quarks except 
in Au+Au collisions at LHC For s quarks, the difference is less than 10% in most collisions, 
which is probably within experimental error and thus not of great interest. 

The PCM reproduces strange particle production in pp collisions at ^/s = 546 GeV 
fairly well, W^ so we assume that it gives reasonable predictions for dN^/dy in AA collisions 



at RHIC and LHC. ||I2[ The factor of 5 increase in s production from RHIC to LHC is 
roughly consistent with most s quarks being produced thermally, which gives dNg/dy oc 
{dN j dyY / A^^^ for AA collisions with the same initial temperature. For S+S collisions 



at SPS, we use data |TB| to estimate dNg/dy. For Au+Au collisions at SPS, we extrapolate 
down from RHIC energy by assuming that dN^/dy oc dN/dy. 



PCM results for c and b production in nuclear collisions [|1^ may be overestimates. 



especially for lower energies, because of the "flavor excitation" processes that are included. 



PCM results for charm production in pp collisions match perturbative results, [14| if only 



flavor creation is included but not flavor excitation. |T^ However, flavor excitation is included 
in PCM simulations of nuclear collisions, so they may overestimate the number of heavy 
quarks that are produced. This overestimate is relatively small at LHC, where the two 
production mechanisms have approximately equal weight, but is large below LHC energy, 
where flavor excitation dominates. 



Using PCM results without flavor excitation from Fig. 4 of Ref. |T2[, dNc/dy\ q ~ 2 for 
central Au+Au collisions at RHIC and 60 at LHC. |16| Simple superposition of perturbative 
pp collision predictions, [|14| assuming that dN^/dy scales as A'^^^ and that the c quarks are 



spread over 4 units of rapidity at RHIC and 5 units of rapidity at LHC, gives respectively 3 
and 10. These rapidity intervals are compatible with PCM results, but smaller than those 



predicted for pp collisions. |jT7| 



Even without flavor excitation, the PCM predicts a large increase in charm production 
(compared to simple superposition of pp collisions) in Au+Au collisions at LHC but not 
at RHIC. The estimated thermal charm production [§ is dNc/dy\ ^ ~ 2 at RHIC and 
10 at LHC, which is too small to account for all of the excess charm at LHC. We make a 
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conservative estimate, taking perturbative predictions |jT^ and scaling dNc/dy as A"^/^, to 
generate the numbers in Table 2. For b quarks, we obtain all rapidity densities by taking 
perturbative QCD predictions for pp collisions Jl^ and assuming that dNb/dy scales as A"^^^ 
(simple superposition) 



C. Predictions for QQ meson production 

Thermal recombination models have been used to analyze strange particle production 
data from ultrarelativistic nuclear collisions, |^,0 so we only briefly discuss mesons containing 
s quarks, and instead concentrate on c and b quark mesons. The only significant difference 
between our model and others in current use [^|-|^ is the inclusion of the pair chemical 
potential, which adds an extra degree of freedom to the models. However, we expect that 
Pgg ^ psPg in nearly all collisions, so the addition of this extra degree of freedom will give 
little change to the predictions of the model, while reducing its predictive power. Thus, 
while the use of figg may be justified theoretically, it is probably not worthwhile including it 
when analyzing collision data. Other models also often include a baryon chemical potential, 
which could be added to our model in the obvious manner. 

We show predictions for /g, the fraction of quarks of flavor Q that form QQ mesons, 
in Table III. We include all confirmed mesons and baryons, as well as those unconfirmed 
mesons that have masses less than the highest confirmed meson mass for each flavor. |1^ 



For simplicity, we take /i^ = 0. The uncertainties in fc and fh are very large, as they depend 
very strongly on u and Tg. 

For fixed z/ and T^, both /^ and fi, drop from S+S to Au+Au collisions at SPS. The 
magnitude of the drop in /^ is largest for large T^ and large z/, and could provide a constraint 
on the two quantities as decrease varies from a factor of 1.2 to 4.5. The decrease in fc could 
provide a second, tighter, constraint, as fc is between 3 and 30 times smaller in Au+Au 
collisions than in S+S collisions at SPS. 

The energy dependences of fb and fc are also interesting. For all cases considered here, 
fb decreases monotonically with increasing energy. The behavior of fc is non-monotonic, 
so changes are smaller and more theory-dependent. Consequently, fc will probably provide 
useful constraints on u and T^ only if z/ ~ 10 and Tc — 150 MeV, in which case fc increases 
from SPS to LHC. 

IV. CONCLUSIONS 

In nuclear collisions at SPS, RHIC, and LHC, s and c quarks are most likely to freeze 
out during the period of two-phase coexistence (at the QCD transition temperature, Tc) if 
u is large corresponding to a strong QCD phase transition. For moderate values of u all 
heavy quark species are most likely to freeze out from hadronic matter. The b quarks freeze 
out at the lowest temperatures. Since these estimates depend mainly on the temperature 
dependence of the entropy density, but are independent of the exact shape, we expect similar 
conclusions to hold even if the transition is not first order as long as there is a reasonably 
sharp change in the entropy density. The mean number of collisions per quark is essentially 
flavor- independent, and large enough to justify the use of quasi-equilibrium models in heavy 

11 



ion collisions at and above SPS energy. Quasi-equilibrium models are only marginally self- 
consistent for S+S collisions at SPS, and are not expected to work for pp collisions at any 
attainable energy. 

In the absence of exotic multi-heavy baryons, density ratios for the various c and b 
resonances will constrain the local freezeout densities, Pq^/pq, but not the global densities, 
Pq. Thus, determination of pc and p^ may require measurement of heavy to light resonance 
ratios. The local and global densities are predicted to be almost equal for s quarks in all 
collisions considered, and for c quarks in collisions at LHC The local b quark density is 
much higher than the global density in all collisions considered here. 

Measurements of the fractions of c and b quarks forming cc and bb mesons, respectively 
fc and fb, may help to constrain u and Tc. This appears particularly promising in light of 
the strong dependence of all other quantities on the entropy density ratio and on the critical 
temperature. 
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TABLES 

TABLE L Freezeout conditions for s, c, and b quarks. A "c" after Tf indicates that freezeout 
was due to loss of thermal contact, while "w" denotes freezeout due to emission. 











Tf (MeV) 






Tf (fm 


■/c) 




Tc 


(MeV): 




150 




200 




150 




200 




Q 


A 


dN/dy 


u: 5 


10 


5 


10 


5 


10 


5 


10 


s 


32 


85 


150c 


150c 


198c 


200c 


2.6 


2.1 


2.0 


1.3 




197 


1000 


145w 


150w 


162c 


200c 


18 


16 


13 


12 




197 


2000 


150w 


150w 


176w 


200w 


23 


22 


21 


18 




197 


3500 


150w 


150w 


197w 


200w 


28 


27 


25 


23 


c 


32 


85 


140c 


150c 


176c 


200c 


5.8 


3.8 


2.9 


2.0 




197 


1000 


116w 


150w 


138c 


182c 


36 


30 


21 


18 




197 


2000 


134w 


150w 


149w 


198w 


46 


39 


33 


29 




197 


3500 


150w 


150w 


168w 


200w 


55 


49 


41 


35 


b 


32 


85 


110c 


145c 


138c 


182c 


12 


10 


6.0 


5.2 




197 


1000 


91w 


120w 


108c 


143c 


75 


65 


44 


38 




197 


2000 


104w 


138w 


117w 


154w 


98 


84 


70 


60 




197 


3500 


117w 


150w 


131w 


173w 


121 


103 


87 


75 










^Q 


(fm) 






"/ 






s 


32 


85 


2.9 


2.5 


2.6 


1.9 


1 


1 


1 


1 




197 


1000 


9.3 


9.3 


8.2 


8.5 


5 


4 


4 


3 




197 


2000 


9.3 


9.3 


9.3 


9.3 


8 


7 


6 


5 




197 


3500 


9.3 


9.3 


9.3 


9.3 


11 


11 


9 


8 


c 


32 


85 


3.2 


2.5 


2.4 


2.0 


1 


1 


1 


1 




197 


1000 


9.3 


9.3 


8.0 


8.1 


6 


5 


4 


3 




197 


2000 


9.3 


9.3 


9.3 


9.3 


9 


8 


7 


6 




197 


3500 


9.3 


9.3 


9.3 


9.3 


13 


12 


10 


8 


b 


32 


85 


3.1 


3.2 


2.3 


2.3 


2 


1 


1 


1 




197 


1000 


9.3 


9.3 


7.9 


7.9 


7 


4 


4 


4 




197 


2000 


9.3 


9.3 


9.3 


9.3 


13 


11 


8 


7 




197 


3500 


9.3 


9.3 


9.3 


9.3 


19 


16 


13 


11 



14 



TABLE II. Freezeout densities. 













PQ (10-^ 


fm-3) 






Pgh/PQPQ 




Tc 


(MeV) 






150 


200 




150 




200 




Q 


A 


dN/dy 


dNg/dy 


ly: 5 


10 


5 


10 


5 


10 


5 


10 


s 


32 


85 


5 


4200 


5300 


5400 


8600 


1.1 


1.1 


1.1 


1.2 




197 


1000 


25 


890 


1000 


1200 


1300 


1.0 


1.0 


1.0 


1.0 




197 


2000 


50 


1400 


1500 


1600 


1800 


1.0 


1.0 


1.0 


1.0 




197 


3500 


250 


5800 


6000 


6400 


7200 


1.0 


1.0 


1.0 


1.0 


c 


32 


85 


0.03 


11 


18 


22 


33 


30 


41 


37 


45 




197 


1000 


0.3 


5.5 


6.4 


9.2 


11 


3.5 


3.1 


3.4 


3.0 




197 


2000 


1 


14 


17 


20 


23 


2.0 


1.8 


1.7 


1.6 




197 


3500 


5 


59 


66 


79 


93 


1.2 


1.2 


1.2 


1.1 


b 


32 


85 


0.0003 


0.06 


0.06 


0.11 


0.13 


6700 


5600 


8100 


6800 




197 


1000 


0.003 


0.03 


0.03 


0.05 


0.05 


530 


450 


500 


430 




197 


2000 


0.02 


0.13 


0.16 


0.19 


0.22 


100 


89 


75 


65 




197 


3500 


0.4 


2.2 


2.5 


3.0 


3.5 


7.4 


6.4 


5.6 


5.0 



TABLE III. Predicted cc and bb meson fractions. 







/c 








h 






Tc (MeV): 


150 




200 




150 




200 




A dN/dy 


u: 5 


10 


5 


10 


5 


10 


5 


10 


32 85 


0.066 


0.075 


0.029 


0.021 


0.62 


0.16 


0.36 


0.072 


197 1000 


0.021 


0.002 


0.008 


0.001 


0.50 


0.053 


0.20 


0.016 


197 2000 


0.009 


0.004 


0.004 


0.001 


0.20 


0.013 


0.070 


0.005 


197 3500 


0.009 


0.010 


0.005 


0.002 


0.075 


0.008 


0.027 


0.002 



15 



